A new screening function for Coulomb renormalization 
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We introduce a new screening function which is useful for the few-body Coulomb scattering 
problem in "screening and renormalization" scheme. The new renormalization phase factor of the 
screening function is analytically shown. The Yukawa type of the screening potential has been 
used in several decades, we modify it to make more useful. As a concrete example, we compare 
the proton-proton scattering phase shifts calculated from these potentials. The numerical results 
document that high precision calculations of the renormalization are performed by the new screening 
function. 
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When one considers the system of few charged parti- 
cles, one faces, as is well known, serious difficulties in the 
calculation of scattering processes !]. This is caused by 
the long range nature of the Coulomb potential in coordi- 
nate space, or, equivalently, its singularity in momentum 
space. 

In order to overcome these difficulties, we would like 
to mention two approaches. One approach is a mod- 
ified time-dependent scattering theoryQ and another 
approach is a "screening and renormalization" method. 
Many investigations have been done in few-body systems 
by the latter approach. For example, Alt et. al. [sj have 
calculated three-body scattering with charged particles 
by the screening and renormalization scheme. Their cal- 
culations worked out successfully. However, the screening 
radius R used in their calculation is about 600 fin. At 
such a large R-value the screened Coulomb potential is 
no longer smooth requiring a careful treatment and in- 
creased computer resource in memory and elapsed time. 
Therefore, it would be desirable to work with a smaller 
R-value. The purpose of this letter is to investigate a new 
screening function different from the one used in Alt et. 
al. which leads to precise results even at small R- values. 

Before we discuss the new screening function, we would 
like to point out the significance of the renormalization. 
If one considers a bound state numerical calculations can 
be performed by the screening method with an appro- 
priately large R-value and the result will be independent 
of the choice of R. However, for a scattering state the 
situation is completely different. The limit of solutions 
achieved with screened Coulomb potentials for increas- 
ing R- values will not agree with the solutions for a pure 
Coulomb potential. The reason lies in the wrong asymp- 
totic boundary conditions going with a screened Coulomb 
potential. A renomalization method, like the one intro- 
duced by Taylor is necessary. 

Alt et. al. used the Yukawa type screened Coulomb 
potential with the screening radius R: 



V R (r) 



,exp[-(r/fl)] 



(1) 



Here in this paper we have proton-proton scattering in 
mind, therefore, e represents the charge of the proton. 
Now we introduce the new screening functions: 



V n R (r) 



,exp[-(r/i?)"] 



(2) 



Note that eq. reduces to the Yukawa type potential for 
n = 1 and is also going to a sharply truncated potential 
for n — * oo, 



lim V n R (r) 



e 
r 



8{R 



(3) 



where 0{r) is the (9-function. The pure and the screened 
Coulomb potentials for n = 1 to 5 at R — 50 fin are 
shown in Fig. 1. This figure reveals that with large n 
there develops a sharper cut-off and the pure Coulomb 
potential is better represented in the inner region. 
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FIG. 1: Pure and screened Coulomb potentials for n = 1 — 5. 

Next, we evaluate the phase shift for pp scattering in 
the state 1 Sq. As a typical NN force we take the Reid 
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Soft Core potential^ in addition to the Coulomb poten- 
tial. There exists the following relation^ among vari- 
ous phase shifts, which can be gained by regarding the 
asymptotic behavior of the wave function: 



6° = 5 1 



+ <j> R + 0(1/ R) 



(4) 



Here 5 C is the phase shift to the strong and pure Coulomb 
potential, S R to the strong and screened Coulomb poten- 
tial, do the standard Coulomb phase shift obtained from 
argT(l + irf) and <pn the renormalization phase. Further, 
i] = e 2 m/2p is Sommerfeld parameter with m the nu- 
cleon mass and p the relative momentum. According to 
Taylor the renormalization phase 4>r(p) is given as 

<Pr(P) = -il I dr 

J(2 P )-i r 

f r exp[-(r/R) n } - 1 , , .„ , 
Jo r 



= -v 



-0(1 /R n ) 
\n(2pR) - 



0(1/R n ), 



(5) 



where 7 is the Euler number (0.5772 ■ • •). The last line 
in Eq. (5) results using the substitution (r/R) n = s/R. 
For n=l this is the result used by Alt et.al. For n — > 00 
one obtains 4>r = —T]hi(2pR), which is the expression 
related to a sharp cut-off. Note that the definition of 
4>r(p) suffer from the uncertainty of 0(l/i?) in eq. J3J 

The phase shifts 5 C , 5 R and 4>r are shown in Fig|2] 
In the evaluation of the phase shifts one has to study 
the dependencies of the cm. energy E cm , the screening 
radius R and the power n in our new screening function. 
As a measure for the quality of the new screening function 
we introduce 




FIG. 2: The comparison of the phase shifts 8 C , S R and 4>R 
for R = 50 fm. 



|A<5| = I 5 C - S R - r)[j/n - ln(2pR)] 



(6) 



This quantity is plotted in Figs EE In Fig El and FigH 

we show the A5|-dependence of E cm for R = 50 fm and 



500 fm, respectively. \A8\ decreases strongly with E cm . 
These figures clearly show that the calculations for n > 2 
are more precise than for n=l. In view of Fig. 1 we have 
to conjecture that this is due to the better approach of 
the pure Coulomb potential by the screened Coulomb 
potential in the inner region if n is larger than 1. 

In Figs. 5 - 6, we illustrate the R-dependence of \A8\ 
for two fixed -E cm -values of 10 and 100 MeV. Trivially, 
|A<$| decreases with increasing R but again the error is 
significantly smaller if n is larger than 1 in comparison 
to n=l. 

Finally Figs. 7 and 8 give the n-dependence of |A<J| 
for fixed E cm and various R's at E cm = 10 and 100 MeV, 
respectively. It is seen that R = 50 fm and n=2 is a good 
enough choice to perform the calculations. 




FIG. 3: |A5| for R — 50 fm against E cm for various n-values. 




FIG. 4: The same as in Fig. 3 for R = 500 fm. 

In summary, we generalized the Yukawa type screening 
potential adopted by Alt et al. ; 3| to the form given in 
Eq. (2). Already for n=2 the screened Coulomb poten- 
tial is closer to the pure Coulomb potential in the inner 
region. The new expression for the renormalization phase 
4>r is just a bridge between the Yukawa type and sharp 
cut-off potential. The numerical results document that 
high precision calculations can be performed by the new 
screening function for n > 2 and a screening radius as 
small as R=50 fm. 
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FIG. 5: |A<5| for E cm = 10 MeV against R for various n. 




acknowledges the supports of the Theory Group of Re- 
search center for Nuclear Physics in Osaka University. 



t 
t 

[1] 



[2] 
[3] 
[4] 
[5] 



Electronic address: yamagu@rcnp.osaka-u.ac.jp 
Electronic address: kamada@mns.kyutech.ac.jp 
Electronic address: koike@i. hosei.ac.jp 
E. O. Alt and W. Sandhas, in Coulomb Interactions in 
Nuclear and Atomic Few-Body Collisions, edited by F. S. 
Levin and D. A. Micha (Plenum Press, New York, 1996). 
J. D. Dollard, Rocky Mount. J. Math. 1, 5 (1971) ;J. Math. 
Phys. 5, 729 (1964). 

E. O. Alt, A. M. Mukhamedzhanov, M. M. Nishonov and 

A. I. Sattarov Phys. Rev. C65, 064613, (2002). 

J. R. Taylor, Nuovo Cimento 23B, 313, (1974) ;M. D. 

Semon and J. R. Taylor, ibid. 26A, 48, (1975). 

R. V. Reid, Ann. Phys. 50, 411, (1968). 



R(fm) 



FIG. 6: The same as in Fig. 5 for E cm = 100 MeV. 
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FIG. 7: |A<5| for E c 
values. 
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FIG. 8: The same as in Fig. 7 for E cm = 100 MeV. 



